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SYSTEM AND METHOD FOR OPTIMIZING OPEN SHORTEST PATH FIRST 
AGGREGATES AND AUTONOMOUS NETWORK DOMAIN INCORPORATING THE SAME 

TECHNICAL FIELD OF THE INVENTION 

[0001] The present invention is directed, in general, to network 
routing and, more specifically, to a system and method for 
optimizing open shortest path first (OSPF) aggregates and an 
autonomous network domain incorporating the system or the method. 

BACKGROUND OF THE INVENTION 

[0002] OSPF is a widely used protocol for routing within an 
autonomous system domain. OSPF implements a two-level hierarchical 
routing scheme through the deployment of OSPF areas. Each OSPF 
area consists of a collection of subnets interconnected by routers. 
Information about links and subnets within an area is flooded 
throughout the area. As a result, every router connected to an 
area knows the exact topology of the area — this includes the 
subnets and the links connecting routers in the area. However, 
details of an area's topology are not advertized beyond the area's 
borders and are thus hidden from other areas. Instead, subnet 
addresses within each area are aggregated and only aggregates are 
flooded into the rest of the network (thus, making an area's 
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subnets reachable from the remainder of the network) . This task of 
advertizing aggregate information about subnets in an area is 
carried out by area border routers (ABRs) , that is, routers 
attached to two or more areas. 

[0003] OSPF areas and address aggregation enable OSPF to scale 
when networks contain hundreds or thousands of subnets. 
Specifically, they help to optimize router and network resource 
consumption, as follows: 

[0004] 1. Router memory. Routing tables only need to contain 
entries corresponding to aggregate addresses for areas not 
connected to the router. Only for areas directly connected to the 
router, are individual subnet addresses in the area stored in its 
routing table. This leads to smaller routing table sizes and thus 
lower memory requirements at routers. 

[0005] 2. Router processing cycles. The link-state database is 
much smaller since it only includes summary information for subnets 
belonging to areas that are not directly connected to the router. 
Consequently, the cost of the shortest-path calculation decreases 
substantially. 

[0006] 3. Network bandwidth. For subnets within each area, 
only aggregate address information (and not individual subnet 
addresses) is flooded into the rest of the network. As a result, 
the volume of OSPF flooding traffic necessary to synchronize the 
link-state databases is significantly reduced. 
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[0007] However, despite the above-mentioned benefits of address 
aggregation, trade-offs are involved. This is because address 
aggregation necessarily results in loss of information, which, in 
turn, can lead to suboptimal routing. In order to see this, one 
must examine in more depth how routing works in the presence of 
aggregation. 

[0008] During routing, each ABR attaches a weight to each 
aggregate advertized by it. This weight is foundational for 
determining the path used by a router external to the area to reach 
subnets covered by the aggregate. From among all the ABRs 
advertizing the aggregate (with possibly different weights), the 
external router chooses the one (say b) for which the sum of the 
following two is minimum: (1) the length of the shortest path from 
the external router to the border router b, and (2) the weight 
advertized by jb for the aggregate. IP packets from the external 
router to every subnet covered by the aggregate are thus forwarded 
along the shortest path from the external router to b and 
subsequently along the shortest path from jb to the subnet. 
However, for a certain subnet, this path may be suboptimal since 
there may be a shorter path from the external router to the subnet 
through an ABR different from jb. The following example, made with 
reference to FIGURE 1, illustrates this problem. 

[0009] FIGURE 1 illustrates an exemplary autonomous network, 
generally designated 100. The network 100 is illustrated as having 
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areas 0.0.0.0, 0.0.0.1, 0.0.0.2 and 0.0.0.3. The network 100 
includes various routers 111, 112, 113, 114, 115, 116, 117, 118, 
119 interconnecting subnets 121, 122, 123, 124, 125, 126, 127, 128. 
FIGURE 1 also illustrates various subnet addresses and weights of 
each link connecting a pair of routers (e.g., routers 113, 115). 
The ABR b x 113 belongs to the area 0.0.0.1, the ABR b 2 114 to te 
area 0.0.0.2, and the ABRs b 3 115 and b 4 116 to the area 0.0.0.3. 
Subnet addresses in the area 0.0.0.3 can be aggregated to different 
degrees. For instance, the aggregate 10.1.0.0/21 covers all the 
subnets in the area 0.0.0.3. In contrast, 10.1.4.0/22 covers only 
subnets 10.1.4.0/24, 10.1.5.0/24, 10.1.6.0/24 and 10.1.7.0/24, 
while 10.1.2.0/23 covers only subnets 10.1.2.0/24 and 10.1.3.0/24. 
[0010] Suppose one of the aggregate addresses advertized by the 
ABRs of the area 0.0.0.3 is 10.1.4.0/22. Suppose further that each 
ABR assigns a weight to the aggregate which is essentially the 
distance of the furthest component in the aggregate from the 
router. Thus, router b 3 115 advertizes 10.1.4.0/22 with a weight 
of 1100 (the distance of the subnet 10.1.6.0/24 from ABR b 3 115), 
while the router b 4 116 advertizes 10.1.4.0/22 with a weight of 
1250 (the distance of the subnet 10.1.4.0/24 from the ABR b 4 116). 
Thus, the ABR b ± 114 belonging to the area 0.0.0.1 forwards all 
packets to the subnets in 10.1.4.0/22 through the ABR b 3 115, since 
the shortest path to the aggregate through the ABR b 3 115 has a 
length of 100+1100=1200, while the shortest path through the ABR b 4 
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116 has length 200+1250=1450. Note, however, that the path from 
the ABR b 2 114 to the subnets 10.1.6.0/24 and 10.1.7.0/24 passing 
through the ABR b 3 115 has a length 1200 and is suboptimal. The 
shortest path from the ABR b x 114 to the subnets is through the ABR 
b 4 116 whose length is only 400. 

[0011] Further, note that considering different weight 
assignments for the aggregates does not alleviate the problem. The 
root of the problem is that a single ABR is selected by the ABR b 2 
114 for reaching all subnets in 10.1.4.0/22. If the ABR b 4 116 is 
chosen instead of the ABR b 3 115, then the paths from the ABR b 2 114 
to the subnets 10.1.4.0/24 and 10.1.5.0/24 through ABR b 4 116 
become much longer (whose length is 200+1250=1450) than the 
shortest paths to the subnets that pass through the ABR b 3 115 
(whose length is 100+50=150) . 

[0012] From the above example, it follows that address 
aggregation can potentially introduce an error in the optimal route 
calculation which is based on shortest paths. For instance, as 
illustrated in the example, when the aggregate address 10.1.4.0/22 
is advertized with weights 1100 and 1250 by the ABRs b 3 115 and b 4 
116, the router Jbj 111 selects the ABR b 3 114 for the aggregate. 
This results in an error of 800 in the selected path from the ABR 
b 2 114 to each of the subnets 10.1.6.0/24 and 10.1.7.0/24. The 
introduced error is essentially the difference between the length 
of the selected path (1200) and the length of the shortest path 
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from the router b z 111 to each of the subnets (400) . Note that the 
selected paths from the router b 1 111 to the other two subnets 
10.1.4.0/24 and 10.1.5.0/24 are the shortest paths, and so no error 
is introduced due to the subnets 10.1.4.0/24 and 10.1.5.0/24. 
[0013] The primary reason for suboptimal paths being selected 
when subnets are aggregated is that a single weight is used by each 
ABR for all the subnets in the aggregate. A single weight may be 
incapable of accurately capturing the distance of the border router 
to every subnet covered by the aggregate. The problem is further 
exacerbated when the subnets are spread across the area with some 
subnets being closer to one border router while others being closer 
to a different border router. This was precisely the case in 
FIGURE 1, where the aggregate 10.1.4.0/22 spans the subnets 
10.1.5.0/24 (closer to the ABR b 3 115) and 10.1.6.0/24 (closer to 
the ABR b 4 116), the distance between whom is greater than 1000. 
As a result, a single weight advertized for aggregate 10.1.4.0/22 
is not representative of the distance between a border router and 
subnets in 10.1.4.0/22. 

[0014] Accordingly, what is needed in the art is a better way to 
select OSPF aggregates and their associated weights thereby to 
avoid suboptimization of routing in a network. 
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SUMMARY OF THE INVENTION 



[0015] To address the above-discussed deficiencies of the prior 
art, the present invention provides systems and method for 
selecting open shortest path first (OSPF) aggregates and aggregate 
weights for a particular area. In one embodiment, an aggregate 
selecting system includes: (1) a database for containing data 
pertaining to candidate OSPF aggregates and corresponding weights 
and (2) an aggregate selector, associated with the database, that 
selects at least a subset of the OSPF aggregates such that the 
shortest path length between the particular source and destination 
subnets resulting from advertisement of a set of weighted 
aggregates approaches the shortest path length between the 
particular source and destination subnets irrespective of the 
advertisement. In one embodiment, a weight selection system 
includes: (1) a database for containing data pertaining to 
candidate OSPF aggregates and (2) a weight assigner, associated 
with the database, that assigns, for the OSPF aggregates, weights 
based on either an average distance of subnets in the area for a 
particular area border router (ABR) of the area or a search 
heuristic . 

[0016] The present invention therefore introduces a better way 
to optimize transmission of data through a network. In one 

embodiment of the present invention, the aggregate selector selects 
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the at least the subset such that the shortest path length between 
the particular source and destination subnets resulting from 
advertisement of a set of weighted aggregates approaches the 
shortest path length between the particular source and destination 
subnets irrespective of the advertisement. 

[0017] In one embodiment of the present invention, the aggregate 
selector treats errors in the shortest path length as having 
unequal degrees of importance. In a related embodiment, the weight 
assigner treats errors in path lengths in the area as having 
unequal degrees of importance. In alternative embodiments, the 
degree of importance for each error is assumed to be equal. 
[0018] In one embodiment of the present invention, the weight 
assigner employs a search heuristic to assign the weights. In a 
manner to be illustrated and described, the problem of assigning 
weights is shown to be JViP-hard. In such cases, a search heuristic 
is appropriate for determining weights. 

[0019] The foregoing has outlined, rather broadly, preferred and 
alternative features of the present invention so that those skilled 
in the art may better understand the detailed description of the 
invention that follows. Additional features of the invention will 
be described hereinafter that form the subject of the claims of the 
invention. Those skilled in the art should appreciate that they 
can readily use the disclosed conception and specific embodiment as 
a basis for designing or modifying other structures for carrying 
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out the same purposes of the present invention. Those skilled in 
the art should also realize that such equivalent constructions do 
not depart from the spirit and scope of the invention in its 
broadest form. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



[0020] For a more complete understanding of the present 
invention, reference is now made to the following descriptions 
taken in conjunction with the accompanying drawings, in which: 
[0021] FIGURE 1 illustrates an exemplary autonomous network; 
[0022] FIGURE 2 illustrates an aggregate tree for eligible 
aggregates covering subnets in area 0.0.0.3; 

[0023] FIGURE 3 illustrates a procedure ComputeMinError () for 
computing aggregates that minimize tree error; 

[0024] FIGURE 4 illustrates a procedure CombineMinError () for 
combining aggregates for a set of trees; 

[0025] FIGURE 5 illustrates an exemplary autonomous network; 
[0026] FIGURE 6 illustrates a diagram of values of Isp and e for 
each literal; 

[0027] FIGURES 7A and 7B illustrate values of lsp and e for each 
conjunct; 

[0028] FIGURE 8 illustrates a procedure ComputeWeightsCumulative 
() for computing weights; 

[0029] FIGURE 9 illustrates a procedure ComputeWeightsMax () for 
computing a W that satisfies Q; and 

[0030] FIGURE 10 illustrates a procedure ComputeWeightsTwoABR () 
for computing weights for two ABRs . 
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DETAILED DESCRIPTION 
[0031] Referring initially to 

[0032] The primary reason for suboptimal paths being selected 
when subnets are aggregated is that a single weight is used by each 
ABR for all the subnets in the aggregate. A single weight may be 
incapable of accurately capturing the distance of the border router 
to every subnet covered by the aggregate. The problem is further 
exacerbated when the subnets are spread across the area with some 
subnets being closer to one border router while others being closer 
to a different border router. This was precisely the case in the 
aforementioned example, where aggregate 10.1.4.0/22 spans subnets 
10.1.5.0/24 (closer to b 3 ) and 10.1.6.0/24 (closer to b 4 ) , the 
distance between whom is greater than 1000. As a result, a single 
weight advertized for aggregate 10.1.4.0/22 is not representative 
of the distance between a border router and subnets in 10.1.4.0/22. 
[0033] One possible way to reduce the error due to the selection 
of suboptimal paths in the presence of aggregation is to avoid 
aggregating distant subnets that are close to multiple border 
routers. Thus, in the example set forth above in the Background of 
the Invention section, instead of advertizing a single aggregate 
10.1.4.0/22, one can choose to advertize two aggregates 10.1.4.0/23 
(with weights 50 and 1250 at b 3 and b 4f respectively) and 
10.1.6.0/23 (with weights 1100 and 200 at b 3 and b 49 respectively) . 
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This reduces the error in the selected paths to 0. Thus, there is 
a trade-off between the number of aggregates advertized (and thus, 
the size of the routing table) and the error in the selected paths. 
Note that the aggregates advertized by border routers do not have 
to be disjoint. That is, for a pair of advertized aggregates, one 
may contain the other. In such a scenario, the longest match 
property of IP routing causes the more specific aggregate to take 
precedence for route computation to subnets within the aggregate. 
The following example illustrates how, by carefully selecting the 
aggregates as well as the associated weights, all the subnets in 
area 0.0.0.3 in the above example can be advertized using only two 
aggregates and incurring zero error. 

[0034] Consider the autonomous network from FIGURE 1. One way 
to ensure that the error in the selected paths to area 0.0.0.3 j s 
subnets is zero is to have the ABRs advertize the following three 
aggregates: 10.1.4.0/23, 10.1.2.0/23 and 10.1.6.0/23. The reason 
for this is that it is possible to choose weights for each 
aggregate at each border router such that the weight equals the 
distance of the border router to every subnet covered by the 
aggregate. For instance, for aggregate 10.1.2.0/23, weights of 
1050 and 250 at border routers b 3 and b 4r respectively, reflect the 
exact distances of the border routers to subnets in it. 
[0035] However, achieving zero error with only two aggregates is 
more challenging. Note that all the subnets in area 0.0.0.3 can be 
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covered by the single aggregate 10.1.0.0/21 however , this 

aggregate by itself, cannot result in zero error. Another 
possibility is to consider the two disjoint aggregates 10.1.4.0/22 
and 10.1.2.0/23 which cover all the subnets. However, as 
illustrated above, since subnets 10.1.4.0/23 and 10.1.6.0/23 
covered by 10.1.4.0/22 are closer to different routers (and distant 
from each other), this cannot result in zero error either. Thus, 
the key to optimizing the error is to bundle 10.1.4.0/23 into one 
aggregate, and 10.1.2.0/23 and 10.1.6.0/23 into the other. It 
turns out that this can be achieved by advertizing the following 
two aggregates: 10.1.0.0/21 and 10.1.4.0/23. The longest match 
characteristic of IP routing causes the latter aggregate to be used 
for routing to subnets in 10.1.4.0/23 and the former to be used to 
route to subnets in 10.1.2.0/23 and 10.1.6.0/23. 

[0036] One question still remains: what weights must be assigned 
to each aggregate? While this is straightforward for the aggregate 
10.1.4.0/23 (since the two subnets 10.1.4.0/24 and 10.1.5.0/24 
covered by it are at the same distance from any given border 
router), it is somewhat less obvious for the aggregate 10.1.0.0/21. 
Simply assigning the weight to be the distance of the most distant 
subnet may not result in the least error. To see this, suppose b 3 
and b 4 advertize 10.1.0.0/21 with weights 1100 and 1250, 
respectively (the maximum distance to a subnet contained in the 
aggregate) . This causes b 2 to select border router b 3 for subnets 
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in 10. 1.0. 0/21 , and the resulting cumulative error in the selected 
paths from b ± to the subnets is 2* (150-150) +2* (1150-450) +2* (1200- 
400) =3000. On the other hand, a lower cumulative error can be 
achieved if b 3 and b 4 advertize 10.1.0.0/21 with weights 730 and 
570, respectively (the average distance to the subnets contained in 
the aggregate) . In this case, b ± selects border router b 4 to access 
subnets in 10,1.0.0/21, and the cumulative error is 2* (1450- 
150) +2* (450-450) +2* (400-400) =2600, which is lower than the error 
for the earlier mentioned selection of weights. 

[0037] Thus, configuring border routers b 3 and b 4 to advertize 
aggregates 10.1.0.0/21 and 10.1.4.0/23, b 3 with weights 730 and 50, 
and b 4 with weights 570 and 1450, causes the error to reduce to 0. 
This is because b 2 selects ABR b 3 for subnets in 10.1.4.0/23 and ABR 
b 4 for the remaining subnets (that is, subnets in 10.1.0.0/21 but 
not contained in 10.1.4.0/23). Thus, the selected paths after 
aggregation are indeed the shortest paths from b x to the subnets in 
area 0.0.0.3. 

[0038] From the above examples, it follows that the two 
important problems related to address aggregation in OSPF are the 
following: (1) selection of the aggregates to advertize at each 
ABR, and (2) assignment of weights to each advertized aggregate at 
each ABR. These two problems are addressed separately. First, for 
a given k and a given assignment of weights to each aggregate, a 
dynamic programming technique is presented to compute the k best 
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aggregates to advertize such that the cumulative /maximum error in 
the selected paths between source destination subnets is minimized. 
This problem is relevant when there is a limit on the number of 
aggregates that can be advertized in order to bound the routing 
table sizes, number of entries in the link-state database or the 
amount of network traffic due to OSPF advertizements . The 
objective then is to choose the k aggregates to advertize such that 
the selected paths are as close to the shortest paths as possible. 
The technique of the present invention can be extended to solve 
variants in which the minimum number aggregates is desired to be 
computed so that the error in selected paths is less than a certain 
user-specified threshold. 

[0039] Later, the problem of choosing weights for each aggregate 
at each ABR such that the deviation of selected paths from the 
shortest paths is minimized is addressed. Specifically, the 
following problem is solved: for each aggregate, determine an 
assignment of weights to the aggregate at each ABR such that the 
error in the selected paths between source destination subnets is 
minimum. While for certain restricted cases, the above problem can 
be solved in polynomial time, the general problem itself is JV3>- 
hard. Search heuristics are presented for selecting weights for 
both the cumulative error as well as the maximum error minimization 
cases . 

[0040] System Model and Problem Formulation 
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[0041] The network is modeled as an undirected graph. Nodes in 
the graph correspond to either routers or subnets. Edges in the 
graph connect routers with other routers or subnets. A link exists 
between two routers if the two routers can directly exchange IP 
packets without going through an intermediate router (that is, the 
two routers are either connected to the same subnet or are 
connected by a point-to-point link) . A link exists between a 
subnet and a router if the subnet is connected to the router. Each 
link has an associated weight, which is the OSPF weight assigned to 
the link (that is, the interface to which the link is connected 
to) . For simplicity, it is assumed that the link is assigned the 
same weight at both ends — the techniques, however, are applicable 
even if link weights are not symmetric. Table 1 describes the 
notation employed herein. 



Symbol 


Description 


k 


Number of aggregates to be selected 


s ,t 


Source, destination subnets 


R 


Set of OSPF areas 


Ri 


Generic term for OSPF area 


s 


Set of all subnets in autonomous system domain 




Set of subnets in area R i 


B 


Set of all ABRs 


B± 


Set of ABRs for area R ± 


b, c 


Generic letters for an ABR 


A 


Set of aggregates eligible for advertising in 
area R £ 


A, 


Set of aggregates eligible for advertising in 
area i? i 
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x, y 


Generic letters for advertized aggregates 


X, Y 


Generic letters for sets of advertized 

a y y i_ " y ct L-Co 


D a ,t 


Degree of importance of the source destination 
subnet oair 


W x 


Weight assignment function for aggregates xn X 


1st (s,t) 


Length of shortest path between subnets s and 
t 


1st (S f t,X, W x ) 


T ^ r-i rr-h "h r\ -F Qhrirfp^t* n^fh Kpfwppn SllbTlGtS S and 

t when aggregates in X are advertized with 
weights in W x 


Table 1 



[0042] The set of subnets S in the network is partitioned into 
disjoint areas. The set of areas is denoted by R and the set of 
subnets in area R ± eR are denoted by A router is said to be 

attached to area R ± if it is directly connected to a subnet in S ± . 
A router that is attached to two or more areas is called an area 
border router (ABR) . B x is the set of ABRs attached to area R ± . 
In addition to area R x (and possibly other areas), every ABR in B i 
is also attached to a special backbone area. The backbone area 
serves to connect the subnets in the various other areas. lsp(s,t) 
denotes the length of the shortest path between s and t. Here, if 
s and t belong to the same area R ir then the shortest path between 
s and t is defined to be over links in area R ± . In case 5 and t 
belong to distinct areas R ± and R Jf respectively, the shortest path 
between s and t involves two ABRs bGB, and cEB jf and consists of 
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three path segments: the first is the shortest path between s and 
b involving links in R if the second is the shortest path between b 
and c over links in the backbone area and the final segment is the 
shortest path between c and t all of whose links belong to area R 3 . 
Note that Isp can be defined in a similar fashion if either of the 
subnets s and t above are replaced by routers. 

[0043] In OSPF, information relating to links and subnets in an 
area are flooded throughout the area. Consequently, routers 
attached to area R ± have detailed knowledge of i?/ s topology. As 
a result, IP packets originating in any subnet s belonging to area 
R ± destined to a subnet t in the same area are forwarded along the 
shortest path between s and t. However, in order to ensure smaller 
routing table sizes and reduce network traffic overhead, detailed 
information about individual subnets within an area are typically 
not be advertized beyond the area's borders. Instead, area i?/s 
ABRs will typically be configured to advertize a set of aggregates 
X that cover subnets in S x and a separate weight for each aggregate 
in X. W x (x,b) denotes the weight assigned to an aggregate xeX by 
ABR beB 1 . Each ABR in B ± floods in the entire backbone area, every 
aggregate xeX along with the weight assigned to x by it; this 
causes ABRs belonging to every other area to receive x. An ABR 
ceB jf in turn, floods x into area R D with an adjusted weight equal 
to minb€B ± {lsp (b,c) +W x (x,b) } . Thus, a subnet s in S jf in order to 
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reach aggregate x that covers subnets in S, selects a path passing 
through ABR jbeB; for which lsp(s,b)+W x (x,b) is minimum. 
[0044] Note that, due to the longest match property of IP 
routing, the most specific aggregate covering a subnet determines 
the path to the subnet. An aggregate x is said to be more specific 
than an aggregate y if x is contained in y, which xey denotes. 
Thus, for a subnet t in S u if x is the most specific aggregate in 
X that covers t, then a subnet s in S jr in order to reach t, 
selects the path comprising of the shortest path from s to b and 
then from b to t, where beB, is the ABR for which Isp (s ,b) +W x (x,b) 
is minimum. lsp(s,t,X,W x ) denotes the length of this selected path 
from s to t for the set of advertized aggregates X and weight 
assignment W x . Thus, Isp (s ,t ,X ,W X ) =lsp (s ,b) + Isp (b , t) , and the error 
in the selected path is simply Isp (s ,t ,X,W X ) -Isp (s ,t) . When s and 
t belong to the same area, Isp (s ,t ,X,W X ) is defined to be equal to 
lsp(s,t). Note that lsp(s,t,X,W x )=°° if X does not contain an 
aggregate that covers t (the implication here is that t is 
unreachable from s) . 

[0045] Further addressed is the problem of computing the set of 
aggregates X advertized across all the areas and the weight 
assignment function W x such that the error in the selected paths is 
minimized. Note that certain restrictions on X and W x need to be 
imposed to ensure the reachability of remote subnets in a different 
area. First, X should be complete, that is, every subnet in S 
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should be covered by some aggregate in X. The next two 
restrictions serve to ensure that an ABR cannot advertize an 
aggregate covering a subnet in S x unless it belongs to B ± . An 
aggregate x is eligible if all the subnets in covered by it belong 
to a single area. Thus, in the autonomous network of FIGURE 1, 
aggregate 10.1.0.0/21 is eligible, since it only covers subnets in 
area 0.0.0.3; however, aggregate 10.1.0.0/20 is not, since it 
covers subnets 10.1.8.0/24 and 10.1.4.0/24 that belong to different 
areas. Let A denote the set of all eligible aggregates such that 
every aggregate in A covers at least one subnet in S. Note that S 
c A. Further, let Ai c A denote the set of eligible aggregates 
that cover subnets in S+. For purposes of the illustrated 
embodiment, it is required that the set of advertized aggregates X 
c A. Further, it is required that only ABRs in B; advertize 
aggregates in A ± . One way to model this is by requiring that 
W x (x,b)=<*> if x=A ± and jbCBj. 

[0046] Now, two problems to be addressed herein may be set 
forth . 

1. Aggregate selection problem: Given a k and a weight assignment 
function W Af compute a complete set X c A containing at most k 
aggregates such that lsp(s,t,X,W x )-lsp(s,t) is minimum. 

2. Weight selection problem: For an aggregate xeA lf compute a 
weight assignment function W M such that £ S/teSex (2sp(s,t, fxi ,IV /Xjl ) - 
lsp(s,t)) is minimum. 
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[0047] Note that in the above-mentioned two problems, every 
source destination subnet pair is assigned the same degree of 
importance. In other words , in the final error, the error in the 
selected path between every subnet pair is treated equally, that 
is, given an equal degree of importance. However, this is somewhat 
restrictive since minimizing the error in the selected paths for 
certain source destination subnets may be more important. This may 
happen, for instance, for subnet pairs between which there is a 
disproportionately high volume of traffic, or subnet pairs carrying 
high-priority or delay-sensitive traffic like voice. Thus, one can 
consider a xx degree of importance" function D which for a pair of 
subnets s f t returns a real value D(s r t) that reflects the 
importance of minimizing the error in the selected path between 
subnets s and t. Note that D(s r t) can be any arbitrary function of 
the volume/priority of traffic flowing between subnets s and t. 
Subnet pairs for whom the error in the selected path does not 
matter (either due to very low traffic volume or due to low 
priority data traffic) can be assigned low values for D(s,t) or 
even zero. The generalized aggregate and weight selection problems 
incorporating the degrees of importance are then as follows. 
1. Generalized aggregate selection problem: Given a k and a weight 
assignment function W Af compute a complete set X c A containing at 
most k aggregates such that Y,s,tes D ( s *t)* {lsp (s ,t ,X r W A ) -lsp (s ,t)) is 
minimum. 
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2. Generalized weight selection problem: For an aggregate xeA it 
compute a weight assignment function W M such that 
Y, s ,tesD(s,t)* (lsp(s,t,X,WA)-lsp(s,t)) is minimum. 

[0048] In the above problem formulations, the cumulative error 
across all source destination subnet pairs is attempted to be 
minimized. An alternative problem of interest is to minimize the 
maximum error across all the source destination subnets. This 
formulation is similar to the above described problems except that 
instead of minimizing Zs.tesO* the objective is to minimize 
max S/teS () . 

Generalized Aggregate Selection Problem 

[0049] This section presents a dynamic programming technique for 
the generalized aggregate selection problem. The technique of the 
present invention exploits the fact that the containment structure 
of aggregates in A is a set of trees. Defined is the notion of 
error for each tree when certain aggregates in it are chosen, and 
shown is the cumulative error of the selected paths when a set of 
aggregates X is advertized is equal to the sum of the errors of the 
trees when aggregates in X are selected in them. Next presented is 
a dynamic programming technique for selecting k aggregates in a 
single tree that minimize its error. Then, a technique for 
combining the results for the collection of trees to derive the 
final k aggregates that yield the minimum error is set forth. 
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Finally, set forth is how the techniques of the present invention 
can be extended for minimizing the maximum error. 
Aggregate Trees 

[0050] The containment relationship among the eligible 
aggregates in A naturally form a set of trees. For aggregates 
x,yeA r an edge from x to y exists if x covers y and every other 
aggregate z€A that covers y also covers x. FIGURE 2 illustrates an 
aggregate tree for the eligible aggregates that cover subnets in 
area 0.0.0.3 (from FIGURE 1). Observe that since A contains 
aggregates 10.1.4.0/22, 10.1.6.0/23 and 10.1.7.0/24, there is an 
edge from 10.1.6.0/23 to 10.1.7.0/24; however, there is no edge 
from 10.1.4.0/22 to 10.1.7.0/24 in the tree. Also, the internal 
nodes in the aggregate tree have either one or two children, but no 
more than two children. For instance, in FIGURE 2, 10.1.0.0/22 has 
only one child since 10.1.0.0/23 does belong to A (subnet 
10.1.1.0/24 does not exist in the network) . Note that each leaf of 
an aggregate tree is a subnet in S. Further, the root r(T) of tree 
T is basically an aggregate that is not covered by any other 
eligible aggregate. 

[0051] Next defined is the error of a tree T when a set of 
aggregates in it have been selected. Suppose x is an aggregate in 
the tree T ,y is the most specific selected aggregate covering x in 
the tree, and X is the set of selected aggregates. The error of 
the subtree rooted at x is recursively defined as follows: 
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E(x,y,X,W A ) = 



[0052] The first condition applies if x has children and xeX. 
[0053] The second condition applies if x has children and x£X. 
[0054] The third condition applies if x has a leaf and xeX. 
[0055] The fourth condition applies if x has a leaf and x$X. 
[0056] The error for an entire tree T with the set of selected 
aggregates X is then simply E (r ( T) , e,X, W A ) (e denotes the empty 
aggregate that does not cover any other aggregate) . Note that each 
recursive invocation of E on x' s children propagates x as the most 
specific selected aggregate if xeX. Consequently, whenever E is 
invoked for the subtree rooted at an aggregate, x,y is always the 
most specific selected aggregate covering x. As a result, the 
error of a tree is simply the sum of the errors of all the leaf 
subnets in it, where the error of a subnet t is 
Y,ses D ( s f t )*(lsp( s 't,{y},W A )-lsp(s,t)), where y is the most specific 
aggregate in X that covers t. Thus, since every subnet in S is 
contained in one of the trees, the sum of errors of all the trees 
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is essentially the cumulative error in the selected paths, which is 
the metric that should be minimized. 

[0057] Thus, the aggregate selection problem has been reduced to 
the problem of computing a set X containing at most k aggregates 
such that the sum of the errors of all the trees is minimum. This 
is then broken into two subproblems that will now be addressed. A 
dynamic programming technique is first presented to compute for a 
single tree the set of aggregates X containing at most k aggregates 
such that the error of the tree is minimum. Then, the computed 
aggregates that minimize the error of each tree are shown how to be 
combined to yield the set of aggregates the minimize the error of 
the entire collection of trees. 
Computing Optimal Aggregates for a Single Tree 

[0058] A set of recursive equations is now presented for 
computing a tight lower bound on the error of a tree assuming that 
at most k arbitrary aggregates in the tree can be selected. The 
equations form the basis of a dynamic programming technique that 
can be used to compute the k best aggregates to select in order to 
minimize the error of the tree. Suppose that x is an aggregate in 
the tree T and y is the most specific aggregate in the tree 
covering x that has already been selected. Then, the minimum error 
of the subtree rooted at x if at most k aggregates in the subtree 
may be chosen, is as follows. 
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min E(x, y 9 k,W A ) = 



E(x,y,<f>,W A ) ifk = 0 

min{min E(u, y 9 k 9 W A \ min E(u, x 9 k-l 9 W A )} 

if k> 0 and x has a single child u 
minlmin^ {min E(u 9 y, i 9 W A ) + min E ( v > y> k ~ z > W a )\ 
min o</</M ( min E ( u > x > *> W a ) + min 

E{v,xX-\-UW A )} } 
if k > 0 and x has children u, v 
min{E(x, y , (f>, W A \ E(x 9 y 9 {x\ W A )} if k > 0 and x is a leaf 



[0059] The intuition underlying the above set of equations is 
that if &=0, then, since no aggregates in the subtree can be 
selected, the minimum error is simply the error of the subtree when 
no aggregates in it are chosen. In case k>0 and x has children, 
and if X is the set of aggregates in the subtree rooted at x that 
if selected result in the minimum error, then the following hold 
for X: (1) either xEX or x$X, and (2) of the remaining aggregates 
in X, i are in the subtree rooted at its left child and the 
remaining k-i or k-i-l aggregates (depending on whether xeX) are in 
the subtree rooted at its right child. Thus, since the error of 
the subtree with x as root is simply the sum of the errors of its 
left and right subtrees, the minimum error can be computed for (the 
subtree rooted at) x by first computing the minimum error of its 
left and right subtrees for 0<i<k selected aggregates and for the 
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cases when x is either selected or not selected, and then choosing 
the combination with the smallest error. Finally, if £:>0 and x is 
a leaf, then there are only two possible alternatives for selecting 
aggregates in x' s subtree: either to select x or not to select x. 
The minimum error for these two cases then yields the desired 
minimum error. 

[0060] Theorem 1: minE (x r y ,k ,W A ) is equal to the minimum of 
E(x,y,k,W A ) f where X is any arbitrary set containing at most k 
aggregates in the subtree rooted at x. 

[0061] From Theorem 1, it follows that minE (r (T) , e, k, W A ) returns 
the minimum possible error for a tree T when at most k aggregates 
in the tree can toe selected. Procedure ComputeMinError in FIGURE 
3 uses dynamic programming to compute the k aggregates that result 
in the minimum possible error for the subtree rooted at x and y is 
the most specific aggregate covering x that has already been 
selected. The procedure is invoked with arguments that include the 
root aggregate of the tree r (T) , e and k. The key ideas are 
similar to those described earlier for the computation of minE, the 
minimum error for the tree. For instance, if an aggregate x has 
children, then procedure ComputeMinError {x,y,l) recursively 
invokes itself for each of its children for the cases when x is 
selected and when x is not selected. Furthermore, in the case that 
x has two children, the procedure is invoked for each child for all 
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the possibilities for the number of aggregates in each child 
subtree . 

[0062] The only difference is that in addition to the minimum 
error, the procedure also computes the aggregates that are 
responsible for minimizing the tree error. Thus, every invocation 
of procedure ComputeMinError , in addition to returning the minimum 
error for the subtree rooted at x, also returns the set of 
aggregates in the subtree that cause the error to be minimum. This 
set is derived by taking the union of the optimal aggregates for 
the subtrees rooted at its children, and adding {x} to it if 
selecting x is required for minimizing the error (steps 11, 20 and 
35) . Note also that in order to improve computational efficiency, 
the optimal aggregates and the minimum error for the subtree rooted 
at x with y as the most specific aggregate and at most 1 selected 
aggregates are stored in subTree [x,y, 1] . aggregates and 
subTree [x,y,l] .error, respectively. The first invocation of 
ComputeMinError {x,y,l) causes the body of the procedure to be 
executed, but subsequent invocations simply return the previously 
computed and stored values. 
Combining the Aggregates for Set of Trees 

[0063] Suppose there are m aggregate trees T 2 , T 2 , T m . 
Further, let T 2 [j] . [error, aggregates] denote the minimum error and 
the set of at most j aggregates in T ± responsible for minimizing 
T/s error. Then, X i [j], [error, aggregates], the minimum error for 
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the set of trees T ± , ... r T z and the j aggregates that minimize 
their cumulative error can be computed using the result of the 
following theorem. 

[0064] Theorem 2: For the set of trees T 2 , T mf 



m \T{j\error ifi = l 

. [i\£rror 1 min ^_ [/}error+ 7J [j - l\erro?\ otherwis 



[0065] Procedure CombineMinError in FIGURE 4 computes in 
X m [k] . [error, aggregates] the minimum cumulative error and the k 
aggregates that minimize the error for the trees T lf T m . After 

computing the error and aggregates for each individual tree in 
steps 1-5, in each iteration of steps 8-14, the XJjj's are 
computed for increasing values of 1 based on the individual tree 
errors and the Xifjj's computed in the previous iteration (as 
stated in Theorem 2). For each X i [j J , the aggregates are computed 
by taking the union of the aggregates for the X^Jl] and TJj-1] 
that result in the minimum error for X i [j]. 
Time and Space Complexity 

[0066] Suppose that d is the maximum depth of an aggregate tree, 
the number of aggregates in A is N and the number of subnets in S 
is n. Note that for 32-bit IP addresses, d < 32 . Then the time 
complexity of the procedure ComputeMinError can be shown to be 
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0 (n 3 +Ndk 2 ) . The reason for this is that lsp(s,t) r the shortest path 
between subnets s and t needs to be computed for all subnet pairs. 
The time complexity of this step is 0(n 3 ). Also, for each subnet 
x and every aggregate y covering it, one can precompute and store 
Y J ses D ( s ^ t ) *(lsp(s,x, {y} ,W A ) -lsp(s,x) ) , thus enabling this 
information to be accessed in constant time. Further, the body of 
ComputeMinError is executed at most once for each combination of x, 
y and 1, For a specific x, there are at most dk different 
possibilities for y and 1 for which the body of the procedure is 
executed. This is because y has to be an ancestor of x in the tree 
and l<k. Each execution of the body of ComputeMinError makes at 
most 2 2 +l recursive calls, and thus, since there are N possible 
aggregates, the total number of times ComputeMinError is invoked is 
0(Ndk 2 ). As a result, the time complexity of procedure 
ComputeMinError is 0 (n 3 +Ndk 2 ) . Further, the space complexity of the 
procedure is 0(n 2 +Ndk), 0(n 2 ) to store the shortest path and error 
information for subnets, and O(Ndk) to store the error and 
aggregate values for each of the Ndk 2 possible combinations of 
values for x, y and 1. 

[0067] It is fairly straightforward to observe that the three 
for loops spanning steps 8-14 of procedure CombineMinError execute 
0(mk 2 ) steps. Thus, the overall time complexity of the procedure 
is 0 (n 3 +Ndk 2 +mk 2 ) , where the first two terms are the time complexity 
of computing the aggregates that minimize the error for the m 
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trees. Note that even though CombineMinError makes independent 
successive invocations to ComputeMinError (r (Ti) ,e ,j) for j=0, . .., 
k, the results computed in subTree during an invocation are shared 
between the invocations. The space complexity of procedure 
CombineMinError is simply 0(mk) to store the X ± and T ± arrays. 
Minimizing Maximum Error 

[0068] Note that instead of minimizing the cumulative error over 
source destination subnet pairs, the techniques set forth herein 
can be adapted to minimize the maximum error over source 
destination pairs. In order to do this, the error of a tree should 
be redefined to be the maximum error of the leaf subnets in it 
(instead of the sum of errors) . Thus, the recursive definition of 
the error of the subtree rooted at x given that y is the most 
specific selected aggregate covering x in the tree and X is the set 
of selected aggregates, is as follows: 



E&y,X 9 W A ) 



ma^ L(sj)*{lsfax, {x\W A )-ls]{s,x)) 



[0069] 



The first condition applies if x has children, and x$X. 



[0070] 



The second condition applies if x has children, and xeX. 



[0071] 



The third condition applies if x is a leaf and x£X. 



[0072] 



The fourth condition applies if x is a leaf and xeX. 
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[0073] Further, the minimum error of the subtree rooted at x if 
at most k aggregates in the subtree can be chosen (given that y is 
the most specific aggregate in the tree covering x and that has 
already been selected), is as follows: 

E{x,y^W A ) ifk = 0 

min{mini?(w, y 9 k, W A \ min E(u, x, & - 1, W A )} 

if k > 0 and x has a singlechildu 

min E(x, yXW A ) = \ mi^min^ {min E{u 9 y, U W A ) + min£(v, y 9 k- i 9 W A )} }, 

min o<,<*-i {min£(w, x, U W A ) + min£(v, x 9 k-l - i 9 W A )}}} 

if k > 0 and x has childrenu, v 
min{E(x, W A \ E{x 9 y 9 [x\ W A )} if k > 0 and x is a leaf 



[0074] Note that, unlike the cumulative error case, where one 
was interested in the distributing the aggregates amongst the 
subtrees of x rooted at children u and v so that the sum of the 
errors of the subtrees was minimized, for the maximum error case, 
one is now interested in minimizing the maximum of the errors of 
the two subtrees (since the error of the subtree rooted at x is the 
maximum of the errors of its two child subtrees) . Thus, the 
following modifications need to be made to procedure 
ComputeMinError to compute the k aggregates that minimize the 
(maximum) error for the tree: (1) replace 
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Y,sesD(s,t)* {lsp(s,x,{y} ,W A )-lsp(s,x)) in steps 5 and 7 by 
max 36S D(s, t;* {lsp(s,x, {y} ,W A ) -Isp (s ,x)) , and (2) replace 
minErrorl+minError2 in steps 26, 27, 33 and 34 with maxfminErrorl, 
minError2g. Similarly, the following simple modification to 
procedure CombineMinError enables it to compute the minimum error 
of a set of trees for the maximum error case: replace X±_ 
Jl] .error+TJj-1] .error in steps 11 and 12 by max {X^ [1] . error , 
T i [j~l ] . error} . 
Weight Selection Problem 

[0075] In the previous section, for a given weight assignment 
function W Af techniques were set forth for computing the optimal 
set of aggregates X for which the error in the selected paths is 
minimized. However, the final error and set of optimal aggregates 
X are very sensitive to the weight that a border router advertizes 
for each aggregate. Thus, the weight assignment problem is 
important for ensuring that selected paths are of high quality, and 
is the subject of this section. 

[0076] Recall that the weight assignment problem is to compute 
a weight assignment function W M for a single aggregate xEA i such 
that the error in the selected paths from all subnets to 
destination subnets covered by x is minimized. The weight 
assignment function N M assigns a weight to x at each ABR beB^ 
Note that the objective is computing the optimal weights for x 
under the assumption that no other aggregates covering subnets in 
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x are concurrently being advertized. Also, since the aggregate x 
is fixed, the subscript M for W is dropped; thus W(b) denotes the 
weight assigned to x by ABR beB i . 

[0077] Intuitively, since W(b) is supposed to represent the 
distance between b and subnets covered by x, two possible logical 
choices for W(b) are the following: 

1 . max tfx 

2. {l/\x\)1 t€K lsp(b f t) . 

[0078] The first choice is simply the maximum distance of a 
subnet in aggregate x from the border router b f while the latter is 
the average distance of subnets in x from b. Note that, since both 
choices are independent of the source subnets (not covered by x) 
and the error to be minimized, as illustrated in the examples 
below, for most cases, neither choice optimizes the objective error 
function. In the following two examples, it is shown that choosing 
W(b) to be max tfx {lsp (b r t) } minimizes neither the cumulative error 
nor the maximum error. 

[0079] Example 2: consider the network in FIGURE 1. Suppose 
one is interested in computing weights for the aggregate 
10.1.0.0/21 that covers all the subnets in area 0.0.0.3. If each 
border router chose the maximum distance to a subnet in 10.1.0.0/21 
as the weight for it, b 3 assigns 10.1.0.0/21 a weight of 1100 
(distance of b 3 from 10.1.6.0/24) and b 4 assigns to 10.1.0.0/21 a 
weight of 1250 (distance between b 4 and 10.1.6.0/24). 
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Consequently, both subnets 10.1.16.0/24 and 10.1.8.0/24 select the 
path through ABR b 3 to access the subnets in 10.1.0.0/21 which has 
a cumulative error of 2*0+2*700+2*800=3000 (for 10.1.16.0/24) and 
2*0+2*900+2*1000=4800 (for 10.1.8.0/24). In contrast, assigning 
weights 1000 and 500 to 10.1.0.0/21 at ABRs b 3 and b 4 , respectively, 
causes the selected paths to be through b 4 which results in a much 
lower cumulative error 2*1300+2*0+2*0=2600 (for 10.1.16.0/24) and 
2*1100+2*0+2*0=2200 (for 10.1.8.0/24). 

[0080] Example 3: consider the network in FIGURE 5. Suppose 
one is interested in computing weights for the aggregate 
10.1.0.0/21 that covers all the subnets in area 0.0.0.3. If each 
border router chose the maximum distance to a subnet in 10.1.0.0/21 
as the weight for it, b 3 assigns 10.1.0.0/21 a weight of 900 
(distance of b 3 from 10.1.6.0/24) and b 4 assigns to 10.1.0.0/21 a 
weight of 1100 (distance between b 4 and 10.1.4.0/24). 
Consequently, both subnets 10.1.16.0/24 and 10.1.8.0/24 select the 
path through ABR b 3 to access the subnets in 10.1.0.0/21 which has 
a maximum error of max{ 0, 700 }=700 (for 10.1.16.0/24) and 
max{ 0, 900 }=900 (for 10.1.8.0/24). In contrast, assigning weights 
1000 and 500 to 10.1.0.0/21 at ABRs b 3 and b 4 , respectively, causes 
the selected paths to be through b 4 which results in a lower value 
for maximum error max{ 700, 0 }=700 (for 10.1.16.0/24) and 
max{500, 0}=500 (f or 10 . 1 . 8 . 0/24) . 
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[0081] Choosing W (b) = ( 1/ I x | ) £ tfX {lsp (b, t) yields somewhat better 
results, because intuitively this is more representative of the 
distance between b and subnets in x than max t0i f lsp (b, t) } . As a 
matter of fact, setting W(b) to be the average distance of b to 
subnets in x can be shown to minimize the cumulative error for the 
weight selection problem. However, it does not optimize the 
maximum error, as Example 4 illustrates. 

[0082] Example 4: consider the network in FIGURE 1. Suppose 
one is interested in computing weights for the aggregate 
10.1.0.0/21 that covers all the subnets in area 0.0.0.3. If each 
border router chose the average distance to a subnet in 10.1.0.0/21 
as the weight for it, b 3 assigns 10,1.0.0/21 a weight of 730 and b 4 
assigns to 10.1.0.0/21 a weight of 570. Consequently, both subnets 
10.1.16.0/24 and 10.1.8.0/24 select the path through ABR b 4 to 
access the subnets in 10.1.0.0/21 which has a maximum error of 1300 
for 10.1.16.0/24 and 1100 for 10.1.8.0/24. In contrast, assigning 
weights 500 and 1000 to 10.1.0.0/21 at ABRs jb 3 and b 4f respectively, 
causes the selected paths to be through b 3 which results in a lower 
value for maximum error 700 for 10.1.16.0/24 and 900 for 
10.1.8.0/24. 

[0083] The following subsections first show that selecting 
W(b) = (1/ I x| )Y jt ex ls P (bf t) results in the minimum cumulative error and 
is a solution to the weight selection problem. However, the 
generalized weight selection problem that involves minimizing the 
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product of the cumulative error of selected paths and their degrees 
of importance is an JVfP-hard problem. Consequently, search-based 
heuristics are presented to solve the generalized weight selection 
problem and a pseudo-polynomial time technique to solve the weight 
selection problem when the objective is to minimize the maximum 
error. Finally, efficient techniques are shown for the generalized 
weight selection problem when B i contains only two ABRs . 
Problem Formulation 

[0084] Some of the notation can now be simplified and new 
terminology can be introduced to address the weight selection 
problem which is: for an aggregate xeA if compute a weight assignment 
function W r such that Y,ses,te*lsp (s ,t, {x} ,W)-lsp (s,t) is minimum. For 
each source s, the selected paths to subnets covered by x is 
through the ABR bEB 1 for which lsp (s,b) + W(b) is minimum (among all 
the ABRs). The ABR selected for source s is designated by B(s,W). 
Note that for tex, Isp (s,t, {x} ,W)=lsp (s ,B (s,W ) ) +lsp (B (s,W) ,t) . 
Further, suppose e(s r b) denotes the error in the selected paths to 
subnets in x if ABR £ is selected for source s. Thus, 
e(s / b)=Y J teKl^P(^^b)+lsp(b,t)-lsp(s f t). Then, 
e (s r B (s,W) )=£ tex lsp (s,t, {x} f W) -lsp (s,t) , and thus the weight 
selection problem becomes that of computing a weight assignment W 
such that Y,ses e ( S ' B ( s ' w ) ) is minimum. 

[0085] The above problem formulation is for minimizing the 
cumulative error. If one wishes to minimize the maximum error, 
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then e (s ,b)=max t€x { Isp (s,b) +lsp (b, t) -Isp (s, t) } and the weight 
assignment N must be such that max seS e (s,B (s, W) ) is minimum. 
Weight Selection Problem (Cumulative Error) 

[0086] For the cumulative error case, it can be shown that 
choosing W(b) to be the average distance of b to subnets in x 
minimizes the cumulative error in the selected paths between 
sources and destination subnets in x. 

[0087] Theorem 3: the weight assignment function W which 
assigns a weight W(b)=(l/|x|) X toc lsp (b, t) to ABR b results in the 
minimum value for £ S(?5 e (s,B (s, W) ) . 

Generalized Weight Selection Problem (Cumulative Error) 

[0088] For the cumulative error case, 

e (s ,b)=Y,tixlsp ( s ' b ) + ls P (b't) ~ ls P ( s 't) is closely related to the 
criterion for selecting an ABR b for s which is that Isp (s,b) +W(b) 
is minimum (note that Y,tex l3 P ( s r t ) is a constant). However, for the 
generalized cumulative error case, 

e (s,b)=Y,t€x D ( s ' t )* {Isp (s,b)+lsp(b,t)-lsp(s,t) ) and thus e(s f b) can 
be any arbitrary value based on the value of D(s,t). This fact 
that e(s f b) can be any arbitrary value makes the problem of 
computing a weight assignment function W that minimizes 
Y, s es e ( s fB (SfW) ) intractable . 

[0089] Theorem 4: for arbitrary values of e(s,jb; and constant 
E, determining if there exists a weight assignment function W for 

which Y, seS e(s,B(s,W) is Jl^-hard. 
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[0090] FIGURE 8 illustrates a randomized technique for computing 
the weights that minimize the cumulative error. The procedure 
keeps track of the best weight assignment found so far in W min/ and 
in each iteration, assigns random weights between 0 and L to some 
subset of ABRs. Here L is chosen to be a fairly large constant — 
a reasonable value for L is one or two times the maximum distance 
between pairs of subnets. 

Generalized Weight Selection Problem (Maximum Error) 
[0091] Recall that, if one is interested in minimizing the 
maximum error, then e (s,b) =max tfx { Isp (s,b) +lsp (b, t) -Isp (s, t) } and 
the weight assignment W must be such that max seS e (s,B (s, W )) is 
minimum. Thus, one can employ an technique similar to the 
randomized search technique to compute a weight assignment function 
that minimizes the maximum error. The only difference is that 
Y,s6se(S' B (S'W)) should be replaced with max sfS e (s,B(s,W) ) in the 
procedure ComputeWeightsCumulative ( ) . 

[0092] However, one can devise a more efficient pseudo- 
polynomial technique for computing the weight assignment that 
minimizes the maximum error. Suppose there exists a procedure P 
that computes a weight assignment W (if one exists) such that 
max sfS e (s,B (s, W) ) <E for some constant E. Then, a simple procedure 
for computing the weight assignment that minimizes the maximum 
error is as follows: (1) Sort the errors e(s,b) between source ABR 
pairs (let E a Er be the errors in order of increasing value), (2) 
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Repeatedly invoke the procedure P for increasing values of i, until 
P returns a weight assignment W for which max sf5 e (s f B (s f W) ) <E i . 
Thus, E t is the smallest value for which a weight assignment exists 
and represents the minimum possible value for the maximum error. 
Further, W is the weight assignment that minimizes the maximum 
error. Note that instead of considering each E ± sequentially, one 
can also use a binary search procedure to compute the minimum value 
for the maximum error more efficiently. 

[0093] The next task is to develop the procedure P that computes 
a weight assignment W (if there exists one) such that 
max s55 e (s,B (s, W) ) <E for some constant E. The problem of computing 
a W such that the maximum error is at most E is equivalent to 
solving a set of inequalities involving the W(b)s as variables. 
For a source s, let ,R denote the set of ABRs beBi for which 
e(s,b)<E; thus, for the remaining ABRs beBi-R, e (s,b) >E. 
Consequently, since the error for each source can be at most E, the 
computed W must be such that one of the ABRs in R is selected for 
s. For this, W is required to satisfy the following set of 
inequalities: min beR { W(b) +lsp (s,b) } <W(c) +lsp (s, c) for all ceB i -R. 
[0094] Thus, for each source s, the set of inequalities 
described above is obtained. Note that the Ws in the equations are 
variables and the Isps are constants. Suppose Q denotes the set of 
inequalities over all the sources. It is straightforward to 
observe that, for a W, the maximum error is E if and only if W is 
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a solution for the set of equations Q. Thus, one needs to focus on 
computing a W that satisfies the inequalities in Q. Observe that, 
if, for a source s, the set R is empty, then a W for which the set 
of inequalities Q is satisfiable does not exist. The reason for 
this is that for the source, inequalities of the form 
min{ } <W(c) +lsp (s,c) are obtainedwhich cannot be satisfied since 
min{}=l. Also, no equations are generated for a source s if R=b ± 
(that is, the error for the source s is at most E irrespective of 
the chosen ABR) . 

[0095] Procedure ComputeWeightsMax in FIGURE 9 is an iterative 
procedure for computing a W that satisfies Q. In each iteration, 
a new weight assignment W new is computed after substituting the 
previous weight assignment W old for the IV' s only on the LHS of each 
inequality in Q (Steps 3 and 4). Note that each inequality in Q' 
has the form C<W (c)+lsp (s,c) , where C and lsp(s,c) are constants 
and W(c) is a variable. Also, lsp max is the maximum value for 
lsp(s,b) for a source ABR pair. 

[0096] ComputeWeightsMax returns a W that is a solution to Q if 
and only if Q is satisfiable. In order to show this, in the 
following lemmas, it is evident that: (1) for any W that satisfies 
Qr W old (b) <W(b) and (2) if Q is satisfiable, then a Sexists that is 
a solution to Q for which Pb&WfbJzi j * ( jJbj j -1 ) 2 ) *lsp max . Thus, 
since W oM does not decrease between successive iterations and the 
procedure terminates only when a W is found or Pb 2 b i W old (b) becomes 
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greater than ( jjfc^j * ( jJbij-1) 2) *lsp max , ComputeWeightsMax computes W 
correctly. 

[0097] Lemma 1: for every ABR beb ir IVj + l old (b) *^j 0 id (b) . 

[0098] Lemma 2: for every weight assignment W that is a 

solution to Q, W(b) *W old {b) . 

[0099] Lemma 3: if Q is satisfiable, then there exists a weight 
assignment W that is a solution to Q and for which Pb 2 b i W{b) < ( ji^j * 
(jjb i j-l)2)*lsp max , 

[0100] Theorem 5: if Q is satisfiable, then procedure 
ComputeWeightsMax returns a weight assignment N that is a solution 
to Q. 

Generalized Weight Selection Problem (2 ABRs) 

[0101] For the special case when b ± contains only 2 ABRs f b x and 
b 2 , it is possible to devise an efficient technique whose time 
complexity is 0{nloqn) where n is the number of sources. Now comes 
time to address the problem of computing a weight assignment W such 
that the cumulative error £ seS e (s,B (s,W) ) is minimized; however, it 
is straightforward to modify the technique of the present invention 
to minimize the maximum error. 

[0102] The error for a source s is efs,^) if 

W{b x ) +lsp (s , b x ) <W(b 2 ) +lsp (s , b 2 ) and e(s,jb 2 ) otherwise. For source s, 
suppose that v(s) =Isp (s, b 2 ) -Isp (s, b x ) and (without loss of 
generality) let v(sl) , . . . , v (sn) denote the values in the sorted 
order. Also, let s 0 and s n+1 be two dummy sources with zero error 
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to both ABRs, v(s 0 )=-«> and v(s n+1 )=°°. Then, for 0<j<n, if 
v ( s j ) <W(b 2 ) -W{b 2 ) <v ( s j + 1 ) , the cumulative error is 
Zi>j e (Si^i) + Li<j e ( s i,^2) • Thus, the problem of computing a W reduces 
to that of computing the v(Sj), v(s j+1 ) pair for which the cumulative 
error is the smallest. 

[0103] Procedure ComputeWeightsTwoABR in FIGURE 10 computes in 
variable V opt the upper value of the v(s 3 ), v(s j+1 ) pair with the 
minimum cumulative error (that is stored in E opt ) . Variable E in 
the procedure is used to keep track of the cumulative error for the 
v(Sj), v(s j+1 ) pair that is currently under consideration in the for 
loop. Thus, choosing values for W(b 1 ) and W (b 2 ) such that W (b 2 )H 
W (b 2 ) = V opt yields the desired weight assignment function W that 
minimizes the cumulative error. 

[0104] Although the present invention has been described in 
detail, those skilled in the art should understand that they can 
make various changes, substitutions and alterations herein without 
departing from the spirit and scope of the invention in its 
broadest form. 
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